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Abstract 

The problem of optimal power allocation among the relays in a two-hop decode-and-forward 
cooperative relay network with independent Rayleigh fading channels is considered. It is assumed that 
only the relays that decode the source message correctly contribute in data transmission. Moreover, only 
the knowledge of statistical channel state information is available. A new simple closed-form expression 
for the average symbol error probability is derived. Based on this expression, a new power allocation 
method that minimizes the average symbol error probability and takes into account the constraints on the 
total average power of all the relay nodes and maximum instant power of each relay node is developed. 
The corresponding optimization problem is shown to be a convex problem that can be solved using 
interior point methods. However, an approximate closed-form solution is obtained and shown to be 
practically more appealing due to significant complexity reduction. The accuracy of the approximation 
is discussed. Moreover, the so obtained closed-form solution gives additional insights into the optimal 
power allocation problem. Simulation results confirm the improved performance of the proposed power 
allocation scheme as compared to other schemes. 
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I. Introduction 

Cooperative relay networks enjoy the advantages of the multiple-input multiple-output (MIMO) 
systems such as, for example, high data rate and low probability of outage by exploiting the 
inherent spatial diversity without applying multiple antennas at the nodes. In cooperative relay 
networks, after receiving the source message, relay nodes process and then retransmit it to 
the destination. Different cooperation protocols such as decode- and-forward (DF), amplify-and- 
forward (AF), coded cooperation, and compress-and-forward can be used for processing the 
message at the relay nodes fl], E3- The benefits of cooperative relay networks can be further 
exploited by optimal power allocation among the source and relay nodes. Specifically, based on 
the knowledge of the channel state information (CSI) at the relays and/or destination, the system 
performance can be improved by optimally allocating the available power resources among the 
relays EJ- 0. 

Different power allocation schemes have been proposed in the literature [|6]- ITT31 . These 
schemes differ from each other due to the different considerations on the network topology, 
assumptions on the available CSI, use of different cooperation protocols for relay nodes, and 
use of different performance criteria [1]. Most of the existing power allocation methods require 
the knowledge of instantaneous CSI to enable optimal power distribution [[TJ, [|6)- [0. The 
application of such methods is practically limited due to the significant amount of feedback 
needed for transmitting the estimated channel coefficients and/or the power levels of different 
nodes. This overhead problem becomes even more severe when the rate of change of the channel 
fading coefficients is fast. 

In this paper, we aim at avoiding the overhead problem by considering only statistical CSI 
which is easy to obtain. Recently, some power allocation methods based on statistical CSI have 
been proposed [fTOll — [[T4|. The optimal power allocation problem among multiple AF relay nodes 
that minimizes the total power given a required symbol error probability (SEP) at the destination 
is studied in ifTTft . The problem of optimal power distribution in a three node DF relay network 
which aims at minimizing the average SEP is studied in [[T2|. The authors of [13] study the power 
allocation problem in a multi-relay DF cooperative network in which the relay nodes cooperate 
and each relay coherently combines the signals received from previous relays in addition to the 
signal received from the source to minimize the average SEP. The power allocation problem 
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aiming at minimizing the average SEP in a cooperative network consisting of two DF relay 
nodes in Nakamgi-m fading channel has been studied in lfl4l . All of the aforementioned power 
allocation methods are based on minimizing or bounding the asymptotic approximate average 
SEP which is valid at high signal-to-noise ratios (SNRs) and is not applicable at low and moderate 
SNRs. 

In our initial conference contribution ffl~5l . a power allocation method for multi-relay DF 
cooperative network with Rayleigh fading channels that minimizes the exact average SEP has 
been proposed. However, the assumption of correct decoding in relay nodes used in iTTSll limits 
its practical applicability. For obtaining a more practically suitable power allocation method, 
we consider in this paper the case when relay nodes may not be able to decode the source 
signal correctly. We derive the optimal power allocation in a multi node DF relay network with 
Rayliegh fading channels in which only relays which have decoded the source message correctly 
contribute in the data relaying. More specifically, after receiving the source message, only the 
relays which decode the message correctly retransmit it to the destination. A new exact and 
simple closed-form expression for average SEP is derived. Then a new power allocation strategy 
is developed by minimizing the exact average SEP rather than its high SNR approximation under 
the constraints on the total average power of all relays and maximum powers of individual relays. 
Only the knowledge of the average channel gains, i.e., the knowledge of the variances of the 
channel coefficients, is assumed to be available. We show that the corresponding optimization 
problem is convex and, thus, can be solved using the well established interior point methods. In 
order to find better insights into the power allocation problem, we derive an approximate closed- 
form solution to the problem and discuss the accuracy of the approximation used. We also show 
by simulations that the exact numerical and approximate solutions provide close average SEP 
performance and that the proposed power allocation scheme outperforms other schemes. 

The paper is organized as follows. System model is introduced in Section II. A simple closed- 
form expression for the average SEP is derived in Section III, while power allocation (exact and 
approximate) strategies based on SEP minimization are derived in Section IV. Simulation results 
are given in Section V followed by conclusion. All technical proofs are given in Appendix. This 
paper is reproducible research [[T6l and the software needed to generate the simulation results 
will be made available together with the paper. It can be also requested by the reviewers if 
needed. 
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Fig. 1. System Model: A single source communicates with a single destination through TV relay nodes and a direct path. 



II. System Model 

Consider a wireless relay network with a single source communicating with a single destination 
through N relay nodes as it is shown in Fig. [Q The relays are assumed to be half-duplex, that is, 
the relays either transmit or receive the signal at the same frequency at any given time instant. 
Therefore, every data transmission from the source to the destination occurs in two phases. In the 
first phase, the source node transmits its message to the destination and the relay nodes, while 
in the second phase, relay nodes retransmit the source message to the destination. The channels 
between the source and the relay nodes, between the relay nodes and the destination, and the 
direct path are assumed to be flat Rayligh fading and independent from each other. The source 
and relay nodes use the M-phase shift keying (M-PSK) modulation^ for data transmission where 
M is the size of the constellation. Relay nodes use DF cooperation protocol for processing the 
received signal from the source node. The received signal at the destination and the ith relay 
node in the first phase can be expressed, respectively, as 

Us,d = \/PQh s , d x + n s4 (1) 
y s ,i = y/Poh Sii x + n Sti , i — 1, . . . ,N (2) 

where x is the source message of unit power, p is the transmit power of the source node, h s d 
and h s i denote the channel coefficients between the source and the destination and between the 
source and the ith relay node, respectively, and n s ^ and n s ^ are the complex additive white 
Gaussian noises (AWGNs) in the destination and in the ith relay, respectively. Since the channel 

'Note that other types of modulation can be straightforwardly adopted and PSK modulation is considered only because of 
space limitation. 
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fading is Rayleigh distributed, the channel coefficients are modeled as independent complex 
Gaussian random variables with zero mean and variances m St d and m Sji , respectively, for the 
channels between the source and the destination and between the source and the zth relay node. 
The additive noises are zero mean and have variance N . 

After decoding the received signal from the source ©, only the relay nodes which have 
decoded the source message correctly retransmit it to the destination through orthogonal channels 
using time division multiple access (TDMA) or frequency division multiple access (FDMA). By 
means of an ideal cyclic redundancy check (CRC) code applied to the transmitted information 
from the source, relays can determine whether they have decoded the received signal correctly 
or not 02]. Then the probability that ith relay node can decode the received signal correctly 
conditioned on the instantaneous CSI can be expressed as ifTTl 



7T 

A ™2 



M-l. 



f M apSK If, .12 

OiW = 1 - - / e A d9 (3) 



where Qpsk = sin (tt/M). 

Let 4> = (0i, 02, • • • , 0jv) t denote a vector that indicates whether each relay has decoded the 
source message correctly or not. Specifically, <fi(i) = 1, if ith relay node decodes the source 
message correctly, and <p(i) = 0, otherwise. Here (-) T stands for the transpose. In the rest of 
the paper, we refer to the vector as the vector of decoding state at the relay nodes. Since 
4> consists of only binary values, there are in total 2 N different combinations that the vector 
4> can take. Moreover, there is a one-to-one correspondence between the binary representation 
of decimal numbers k = 0, . . . , 2 N — 1 and different values that the vector can take. For 
example, in a cooperative network with two relays, if the relay enumerated as first decodes 
the source message correctly and the relay enumerated as second decodes it incorrectly, the 
corresponding decoding state vector is (1,0) T and the corresponding representation in decimal 
is 1 ■ 2 1 + ■ 2° = 2. For simplicity, we represent hereafter each possible combination of vector 
4> by its corresponding decimal number k and denote this combination as <fi k . For example, O 
corresponds to the situation when all relay nodes decode the source message incorrectly. 

The received signal from the ith relay node which is able to decode the source message 
correctly, that is, (f>(i) = 1, at the destination can be modeled as 

Vi.d = VPihi, d x + n i4 (4) 
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where Pi is the transmitted power of the ith relay node, h id is the channel coefficient between the 
ith relay node and the destination, and n iyd is the AWGN with zero mean and variance N . The 
channel coefficient h id is modeled as an independent complex Gaussian random variable with 
zero mean and variance m id due to the Rayleigh fading assumption. It is worth stressing that 
the assumption that the channel coefficients are independent from each other is applicable for 
relay networks because the distances between different relay nodes are typically large enough. 
It is assumed that the destination knows perfectly the instantaneous CSI from the relays to the 
destination and the instantaneous CSI of the direct link. The knowledge of instantaneous CSI 
for the links between the source and the relay nodes is not needed. Then the maximal ratio 
combining (MRC) principle can be used at the destination to combine received signals form the 
source and relay nodes. As a result of MRC, the received SNR at the destination conditioned 
on the decoding state at the relay nodes, i.e., <fi k , can be expressed as 



where 7 S = Po\h S)d \ 2 /N and ji = Pi\h itd \ 2 / N are the received SNRs from the source and the zth 
relay node at the destination, respectively. Here 7 S and 7$, i = 1, • • • , N are exponential random 
variables with means pom Sid /N and pim iid /N , respectively. Moreover, 7^ and %, i = 1, ■ ■ ■ ,N 
are all statically independent. 



For the considered case when the data transmission is performed using the M-PSK modulation, 
the SEP of the signal at the destination conditioned on the channel states CS = {h s>d , h s i , h ijd } 
and the decoding state at the relay nodes <p k can be written as ifPTl 



Using the total probability rule, the SEP conditioned on the channel states can be expressed as 





(5) 



{i\<l> k (i)=l} 



III. Average SEP 




(6) 




(7) 



k=0 



where P r {<p k } is the probability of the decoding state <p k that can be calculated as 



n (!-«*) 



(8) 



Wfc(0=i} 



{i|* fc (i)=0} 
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where is the probability of correct decoding in ith relay node ©. Note that for obtaining ©, 
the independency of the AWGNs at the relay nodes has been exploited. 

The average SEP can be obtained by averaging © over h s d , h Sji , i = 1, ■ ■ • ,N, and h i:d , 
i = 1, • ■ ■ , N and using the fact that P r {4> k } is statistically independent from P e {CS, <p k }. The 
latter follows from the statistical independence between the channel coefficients and the fact that 
P r {4> k } depends only on {h Sti }f =1 and P e {CS,<p k } depends only on {h it d}fLi an d h S:d . Then 
the average SEP can be expressed as 



where bi = gpsK^i^/No, i — 1, . . . , N and b = gpsKiTi St d/N . Similarly, the first expectation 
in © can be computed as 




(9) 



where E{-} denotes the expectation operation. 
The second expectation in ([9]) can be computed as 




(10) 




a • n c 1 - ^ 



(id 



{i|# fc (i)=l} {i|* fe (j)=0} 



where 




(12) 



and a = gpsKm Si i/N Q , i = 1, . . . , N. 
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Substituting ( flOl ) and (fTTT) in ©, the average SEP can be equivalently expressed as 



2™-l 



P. 



^e n a n (i-A) 



X 



fc =o \{*l* k (i)=i} OI«A/bW=o} 

sin 2 (0) 



A/ 



n 



sin (0) + 



d0 



(13) 



{i\$ k (i)=l}U{i=0} 

Setting /3 = 1 for notation simplicity, the average SEP expression in (TT3T ) can be simplified as 



ft 



x 




(14) 



e n a-A) 

fc=0 {i|0 fc (i)=O} 



X 



n 



A ■ 



sin 2 (0) 



sin (0) + hp, 



d6 



(15) 



W fe (i)=l}U{i=0} 

Finally, the expression (TT3T) can be further simplified as 



P« = - 



7T 



M^lir N 



i=0 



/W(0) 



sin (0) + hpi 



de. 



(16) 



To verify the latter result, note that (Q~5) can be obtained by simply expanding (fT6l) . Moreover, 
after finding the integral in (fl"2l ). fy, i = 1, . . . , N can be expressed as 



QPo 



tan" 1 (cot (^pvr) 



CiPO+1 



QPo + 1 
1 tan" 1 (cot (^tt)) 



V 



IT 



N. 



(17) 



2 7T 

Here ft has a meaning of statistical average of the correct decoding probability in the ith 
relay node @ with respect to the channel coefficients. It is worth noting that the average SEP 
expression (fT6l) is used later for finding the optimal power distribution among the relay nodes. 

Using partial fraction decomposition, it is possible to further simplify (fT6l) and find a closed- 
form expression for the average SEP. It is worth noting that for simplicity and because of space 
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limitations, we hereafter assume that bi ^ bj, i ^ j, i,j — 0, . . . , N. However, if the condition 
bi 7^ bj, i 7^ j, i,j = 0, . . . , N does not hold, the following average SEP derivation approach 
remains unchanged, while the only change is the need of using another form of the partial fraction 
decomposition. Thus, it is straightforward to derive closed-form expression for the average SEP 
in the case when hi ^ bj, i ^ j, i,j — 0, . . . , N does not hold by using the same steps as we 
show next. 

Toward this end, we first rewrite the integral inside dT4l >. denoted hereafter as /, by changing 
the variable as 

Sin 2 (g) 

{i|</> fc (i)=l}U{*=0} v > Wl 



M ' 







II fl J_h. n . \+h.*n.r2 



ni ( m-i \ 1+x 2 - LJ - (1 + biPi) +biPiX 2 

"'HtT 1 ] {j|0 fe (i)=l}u{j=O} V %yiJ %F% 



(18) 



Then applying the partial fraction technique to the right-hand side of ([TBI and using the fact 
that bi bj, i ^ j, i,j — 0, . . . , N, we obtain 

1 



n 



(b'l*fc(j')=l}U{j'=0})-{j#i} x "jPj 

Y 7 r~, 5 —dx. 



1 + — ) + X 

{i|<£ fe (i) = l}U{i=0} V b iP 



(19) 



The integral in (TT~9T > is the summation of the terms that all have the form of l/(a+ x 2 ), a > 0. 
By using the fact that J l/(a + x 2 )cb = l/v/titan" 1 ^/^), a > 0, (fl9l) can be calculated as 



fi|<Mi)=l}U{i=0}> 



tan COt r^7T 



M-l \ / \ 7T 



M JVb iPi +l 2 



X 



F (O1^ fc (i)=i}u{j=o})-{i^} 1 6 iPi 



+ ( - - tan" 1 (cot (^—^tt)) ) (20) 



2 M 

Substituting the so obtained expression (l20l) in (IT41) and also expanding the resulted expression 
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the average SEP can be equivalently expressed as 

2^-1 N 

p *= e n ft- n ( x -ft) 

A:=0 {i|0 fe (i)=l}U{i=O} {i|<^fc(i)=0} 



{i|</. fc (i)=i}u{i=o} 



W^)Vb^.) 1 



7T 2 



n 



biPi + l J-- 1 - i_£z£l 

({7'l^(i)=i}u&=o})-{,Vi} 



2-^-1 N 



+ e n ax n ^-ft) 

fc=0{i|</, fc (i)=l}U{i=0} {#fc(*)=0} 

1 tan" 1 (cot (^-tt)) 



(21) 



2 7T 

Moreover, we can modify the first term of (I2T1 by taking the product ri{i|^ fc (i)=i}u{i=o} ft mto 
the summation part of that term and expanding the product such that is multiplied by the 
corresponding term in the summation and also modify the second term of (|2T|) by using the 
following equality which is easy to prove 

2 N -1 N 

e n ft - n (i-A)=i- (22) 

fc=0{i|^ fc (i)=l}u{i=0} {i|4> fc (i)=0} 

By doing so, the average SEP can be rewritten as 

2 N -1 

p * = e n ^-ft) 

fc=0 {i|* fe (i)=0} 



x E (ft 

{i|<Mi)=i}u{;=o} 



'tan-Hcot^Tr)^^) ! 

7T 2 



X 



n 



F ({7l*fcO>l}U{7'=0})-{?'Al 1 



+ A tan" 1 (cot(^7r)) \ _ ^ 
Finally, after rearranging and factoring out the terms, we can obtain the following closed form 
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expression for the average SEP 




tan-Hcot^vr)^) ! 



) 



biPi 



" n (Ak^-m 





1 tan" 1 (cot (^r"7r)) 



(24) 



Note that (|23l) can be obtained by simply expanding (1241) . It is interesting to mention that if we 
consider similar conditions as in [fT51 , in which all relays are capable of decoding the source 
message correctly, i.e., $ = 1,2 = 1, . . . , N, and also there is no direct link between source 
and destination, i.e., b = 0, the closed-form average SEP (1241) simplifies to the one that was 
obtained in [15]. The closed-form expression (1241 is simple and does not include any other 
functions rather than basic analytic functions. 



In this section, we address in details the problem of optimal power allocation among the relay 
nodes such that the average SEP is minimized. Only statistical information on the channel states 
is used. Moreover, we assume that the power of the source node po is fixed and, in turn, the 
average probabilities of the correct decoding of the relays, i.e., i = 1, . . . , N, are also fixed. 
The relevant figure of merit for the performance of relay network is then the average SEP that 
is derived above in closed form. It enables us to apply power allocation also in the case when 
the rate of change of channel fading is high. Since only statistical CSI is available and relay 
nodes retransmit only if they decode the source message correctly, the averaged power of the 
relay nodes over CSI and probability of correct decoding need to be considered. It is easy to see 
that the average power used by the ith relay node equals fopi. Indeed, since (3i is the statistical 
average of the probability of correct decoding in the ith relay node, the transmitted power of 
the zth relay node pi weighted by gives average power used by the ith relay node during a 
single transmission. 

In the following, we develop a power allocation strategy by minimizing the average SEP 
(fT6l) . while satisfying the constraint on the total average power used in relay nodes pn and 
the constraints on maximum instantaneous powers per every relay p™ ax , i — 1, . . . , N. With the 
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knowledge of the average channel gains m Sii = E{\h s ^\ 2 }, m^ d = E{\h itd \ 2 }, m Sjd = E{\h Sjd \ 2 }, 
and specifications on p R , p™ ax , i = 1, . . .,N, and p , the power allocation problem can be 
formulated as 

min P e (p) (25) 
p 

N 

s.t. ^2(3 iPi = p R (26) 

0<p i <p?°* : ,i = l,...,N. (27) 

where P e (p) is the average SEP (fT6l) written as a function of powers at the relay nodes p = 
(pi, . . . ,Pn) T ■ For notation simplicity, we use hereafter the following equation for the average 
SEP in CU) 

A/-1 _ 
1 f — 



where 



Pe(p) = ~ 9(0,p)d6 (28) 

p)-nf(i-A)+A -^g\ y (29) 

Note that the optimization problem (|25T>— (|27T> is infeasible if X)i=i P>iVT ax < PR anc ^ ^ has the 
trivial solution, that is, pi = p™ ax , if X)i=i PiPT ax = PR- Thus, it is assumed that J2iLi PiPT ax > 
Pr. The following theorem about convexity of the optimization problem (I25l-(l27l is in order. 

Theorem 1: The optimization problem ri25l)-rf27l) is convex. 

See the proof in Appendix A. 

Although this problem does not have a simple closed-form solution, an accurate approximate 
closed-form solution can be found as it is shown in the rest of the paper. A numerical algorithm 
for finding the exact solution of the optimization problem (l25l) -(l27T) based on the interior-point 
methods is summarized in Appendix B. Despite the higher complexity of the numerical method 
as compared to the approximate closed-form solution, numerical method can provide an exact 
solution, which can be used as a benchmark to evaluate the accuracy of the approximate solution. 



A. Approximate Closed-Form Solution 

The optimization problem (|25l) -(l27l) is strictly feasible because Yl!i=\ PiPT ax > PR as ^ has 
been assumed earlier. Thus, the Slater's condition holds and since the problem is convex, the 
Karush-Kuhn-Tucker (KKT) conditions are the necessary and sufficient optimality conditions 
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lfT8l . Indeed, since Yli=iPiPT ax > PR' men Pi = PT ax Pn/ 'Q_)i=i PiPT ax ) * s a strictly feasible 
point for the optimization problem (|25T)-(|27T). 
Let us introduce the Lagrangian 



1 /'~ 7r \ - 

£(p, \ 7, f ) =- / g(9, P)d9-^ XiPi 
71 J o i= l 

N / N 

+5^7i(Pi-pr a )+H ^te-^l (30) 



i=l 

where A = (Ai, • • • , Aat) t and 7 = (71, • • • , 7jv) t are AT x 1 vectors of non-negative Lagrange 
multipliers associated with the inequality constraints p, L > 0, i = 1, • • ■ , N and pi < p™ ax , i = 
1, • • • ,N, respectively, and v is the Lagrange multiplier associated with the equality constraint 
Yli=i Pi — PR- Then the KKT conditions can be obtained as 

Ai>0, 7i >0, z = l,-.-,iV (31) 

XiPi = 0, vfa - pT ax ) = 0, i = 1, ■ ■ • , N (32) 

1 f^ n -/3 i b i sm 2 (6)g(e,p) 
%J sin 4 ( 



de 



i\6) + (2 - pjbffii sin 2 (#) + (l - /3i)b1pf 
+^A-Ai + 7 = 0, i = l,---,N (33) 

0< Pt <pT ax , i = l,---,N (34) 

N 
i=l 

Although the exact optimal solution of the problem (l2~5T)-(l27T) can be found through solving 
the system (|3D)-(|35T) numerically or as it is summarized in Appendix B by solving the original 
problem directly, a near optimum closed-form solution for the system (|3TI)-(|35|). and thus, 
optimization problem (l25T)-(l27T) can be found by approximating the gradient of the Lagrangian, 
that is, the left hand side of d33l) . Specifically, it can be verified that for fixed p^, i = 1, . . . , N 
the function g(9,p) is strictly increasing/decreasing with respect to 9 in the intervals (0, 7r/2) 
and (tt/2, (M — l)n/M), correspondingly. Under the condition that the number of relays is large 
enough, the slope of the increment and decrement in the aforementioned intervals is high and 
the Chebyshev-type approximation on the conditions (|33l) is highly accurate. This approximation 
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of the conditions ( 1331 is of the form 

PA 



7r(l+(2-/3 i )6^+(l-A)^f) 

M-l 

x/ g(e,p)d0 + friz- Ai + 7i = 0, 
■/ o 

i = l,---,N (36) 

where the fact that sin 2 (7r/2) = 1 is used in the ratio 

-PilH sin 2 (fl) 

(37) 



sin 4 (0) + (2 - A)^sin 2 (^) + (1 - P^pf 



By rearranging the denominator of ([371 ). substituting the approximation (1361 ) in (1331 ). dividing 

M-l 

the equations (13Tb— (f33l) by the positive quantity J Q M p(0, p) d6, and also multiplying these 
equations by ir, we obtain 

a:>o, 7 ;>o, i = i,..-,iv (38) 

Afo = 0, -rtfa - vf ax ) = 0, i = 1, • • ■ , JV (39) 

~A(i+^) + (i-A)(i+^) 2 A -^ +7i " u 

i = l,-..,JV (40) 

0<Pi<pr x , i = 1, - - • , (41) 

JV 

^PiPi = PR- (42) 

8=1 

M-l . M-l . M-l 

where z/ = i/tt/ J q ~ g{6) dO, ^ = ^n/ J ~ g(9) d9, and A^ = \ix / J Q ~ g{9) dO, i = 
1, ■ • ■ , N. It is possible to eliminate A- from the set of equations (l38l)-(l42l in order to find a 
simpler set of smaller number of equations. By doing so, the approximate KKT conditions can 
be equivalently rewritten as 
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7 />0, i = l,...,N (43) 

l/(p l -p l max ) = 0, i = l,---,N (44) 
fa/, i & b i ^ - 



j 9 i (l+6 iPi ) + (l-A)(l+6<Pi)V ' 

z = l,-. -,JV (45) 



■ < p>' + 7/, 



A(l+6«Pi) + (l-A)(l+6iPi) 2 

z = l,-..,JV (46) 

0< Pi <Pi max , i = l,--- ,N (47) 

JV 
i=l 

Note that (|40l) is rewritten as ([461 since Aj or, equivalently, A- is positive. Moreover, (|45l is 
obtained by solving (|40l) with respect to A. and then substituting the result in A^jOj = 0. The 
following result gives a closed-form solution for the system (I43l)-(l48l). and thus, it also gives an 
approximate solution for the power allocation optimization problem (|25l-(|27l. 

Theorem 2: For a set of DF relays, the approximate power allocation {pi , . . . , Pn}> i-e-, 
the approximate solution of the optimization problem rf25l)-rf27l). is 

26,(1-/3.) h ] 



o 



where v' is determined so that J2iLi PiPi = PR- 
See the proof in Appendix C. 

It is interesting to mention that the power allocation scheme of [fT5l is a special case of the 
power allocation method given by (|49l when 0i — > 1. Indeed, it can be checked that for — > 1, 
we have 

-A+yi+4]H^V l 

Pi = hm 



ft-n 26,(1-/3,) 6, 
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(50) 

/ o 

This special case shows that the total power should be distributed only among relays from 
a selected set and the relays with 'better' channel conditions use bigger portion of the total 
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power. Note that this solution can be viewed as a form of water-filling solution and has similar 
complexity. However, the solution (|49l ) is not interpretable only by the means of water-filling, 
and the optimal power allocation among the relays and relay admission depend on the average 
probability of correct decoding by relay nodes and also on the corresponding average channel 
gain-to-noise ratios for admitted relays. 

V. Simulation results 

We consider a cooperative relay network consisting of a source-destination pair and N = 5 DF 
relay nodes. In the first phase, source transmits its message to the destination and the relay nodes, 
while in the second phase, only the relay nodes, which decoded the source message correctly, 
retransmit it to the destination. All the nodes use QPSK modulation for data transmission and 
noise power is assumed to be equal to 1. Relays are assumed to be located on a line of unit 
length connecting the source and the destination nodes. Positions of the relay nodes with respect 
to the source node are then randomly selected according to the uniform distribution and are 
{0.0117, 0.1365, 0.2844, 0.4692, 0.8938}. All the channel coefficients are modeled as zero mean 
complex Gaussian random variables. The variance of the channel coefficient between pth and 
gth nodes is l/c^ where d PA is the distance between the nodes and 2 < v < 6 is the path-loss 
parameter which is assumed to be equal to 3 throughout our simulations. 

A. Closed-Form SEP 

We first aim at comparing the closed form SEP expression ([24l with the SEP obtained through 
Monte-Carlo simulations and demonstrating their equivalence. The source node transmits its 
message to the destination node through only the relay nodes 1, 3, 5 or all relay nodes 1, 2, 3, 4, 5. 
The total power is equally divided among the source and the relay nodes. Fig. 1 shows the average 
SEP corresponding to the closed-form expression (l24l and the SEP found by Monte-Carlo 
simulations with each point obtained by averaging over 10 8 independent runs. The corresponding 
SEPs are shown versus the total power. Fig. 1 confirms the fact that the closed-form expression 
(1241) results in the same average SEP as the one obtained through numerical simulations. The 
closed-form-based SEP and the SEP obtained through numerical simulations coincide in both 
scenarios considered with 3 and 5 relay nodes. 
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B. Accuracy of the Approximate Closed-Form Solution 

In this simulation example, we study the accuracy of the proposed approximate power alloca- 
tion method (see Theorem 2) by comparing its performance to that of the optimal SEP found by 
solving the exact problem (l25T)-(l2"7T) using the algorithm summarized in Appendix B. Moreover, 
the performance of the proposed power allocation strategy is compared with that of the equal 
power allocation based strategy. The source transmit power is assumed to be fixed and equals 
to 1. The power transmitted by each relay node in the scheme with equal power allocation is 
Pi = PrJ Yld=i = 1, ■ ■ ■ , R- For evaluating the SEP, the closed form the SEP expression in 
(|24|) is used. 

Figs. [3] and |4] illustrate the SEP of the aforementioned power allocation methods when only 
the relays 1, 3, 5 or all relays 1, 2, 3, 4, 5 are used for the cases when p™ ax — p r , % — 1, . . . , N and 
pmax _ p r /2, i = 1, . . . , N . From these figures, it can be observed that the SEP corresponding to 
the power allocation obtained by approximating KKT conditions (|49| ) is very close to the optimal 
SEP obtained through numerical solution of the optimization problem (l25T)-(l27T) even in the case 
when there are only 3 relay nodes. It is noteworthy to mention that our extensive simulation 
results confirm the high accuracy of the proposed approximate power allocation method for 
different number of relay nodes. It can be observed from Figs. [3] and @] that there is more than 
5 dB performance improvement for the proposed optimal power allocation scheme as compared 
to the equal power allocation scheme in high SNRs and when, p™ ax = p r . As it is expected, if 
maximum available/allowable power at the relay nodes is limited, the corresponding performance 
improvement deteriorates. 

C Performance Improvement Compared to Recent Methods 

In our last simulation example, we aim at comparing the performance of our proposed power 
allocation scheme with another recently proposed power allocation scheme. The other scheme 
aims at minimizing the power required by all AF relay nodes under the constraint that the average 
SEP does not exceed a certain desired value and it uses the knowledge of the average channel 
gains ifTTf. Thus, as compared to our proposed power allocation scheme, which minimizes the 
average SEP of a DF cooperative network subject to a fixed total power of the relays nodes, 
the method of [fTTI minimizes the total power required by the relay nodes of an AF cooperative 
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network and fix the bound on the average SEP. Both schemes use only the knowledge of the 
average channel gains. 

Compared to the method of lITTI according to which all relays always retransmit the source 
message, our proposed method exploits the additional information about the source message 
decoding failure/success in each relay node, and thus, not all relays retransmit the message. As 
a result, it is expected that the proposed method will outperform the one in 0T). In addition, an 
approximation of SEP, which is accurate only at high SNRs, is used in 0~T|, while the exact SEP 
expression is derived and used for our power allocation scheme. To ensure a fair comparison 
between two aforementioned power allocation schemes, we first find optimal power allocation 
according to the scheme of [|TT|. Then, we use the total power obtained by the scheme of |fTTfl 
as a power bound for our proposed power allocation scheme to ensure that the network operates 
with the same amount of total power in both cases. 

Fig. [5] compares the SEP corresponding to the optimal power allocation of the AF cooperative 
network using [IITI and the SEP corresponding to the optimal power allocation obtained based 
on the proposed method for DF networks for relays 1, 3, 5 and all relays 1, 2, 3, 4, 5, respectively, 
versus different desired SEPs denoted as SEP for the case that there is no restriction on the 
maximum allowed power for each relay node. Monte-Carlo simulations with 10 6 independent 
runs are used for obtaining each point in Fig. |5] It can be observed from this figure that, as 
expected, the proposed power allocation scheme has superior performance over that of the scheme 
of ifTTIl with the same average transmit power from the relay node. This improvement can be 
attributed to the fact that in the proposed scheme the relay for which the channel condition 
between the source and relay is not good enough or, equivalently, the relay that can not decode 
the source message correctly, does not transmit. This additional information is exploited in the 
proposed scheme, while it is not used in the scheme of [11]. In addition, the exact average SEP 
is used for the proposed scheme versus the approximate one in the scheme of 0TJ. 

VI. Conclusion 

A new simple closed form expression for the average SEP for DF cooperative relay network 
has been derived. Using this expression, a new power allocation scheme for DF relay networks 
has been developed by minimizing the average SEP under the constraints on the total power of 
all relays and the maximum powers of individual relays. The proposed scheme requires only the 
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knowledge of the average of the channel gains from relays to destination and the direct link. 
The exact and approximate closed form solutions to the corresponding optimization problem 
have been found and a high accuracy of an approximate solution is demonstrated. According 
to the proposed approximate solution, the power allocation and relay admission depend on the 
average probability of correct decoding by relay nodes and also on the corresponding average 
channel gain-to-noise ratios for admitted relays. The improved performance of the proposed 
scheme compared to some other schemes is convincingly shown via simulations. 

Appendix A 

Let us consider the following formula for the average SEP 
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The Hessian of the integral inside (ISTI) with respect to p can be obtained as 
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D*(*)(0»p) = 

l (sm 2 (0)+6 1 p 1 ) 2 '"''(Bin 2 (0)+6^)V 
with denoting the Hermitian transpose and diag(-) standing for a diagonal matrix. 

Since the matrices a^ fc (^, p)a^ (0, p) and D<£ fc (6>,p) are both positive semi-definite for 9 £ 
(0, (M — 1)tt/M) and p» > 0, z = 1, • • • , iV and also 5V> fc (0,p) > 0, the Hessian matrix is 
positive semi-definite. The average SEP (15TT) is a linear combination of integral expressions that 
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are convex, and as a result, it is convex on nonnegative orthant. Moreover, since the constraints 
of (T25T)-(l27l) are linear and form a convex set which is a subset of nonnegative orthant, the 
problem (l25l)-(|2"7T) overall is convex. □ 

Appendix B 

The numerical procedure for solving the problem d25li-(l27Ti is based on the interior-point 
methods. Specifically, the barrier function method, which is one of the widely used interior- 
point methods, is applied. The barrier function method for solving the problem (l25l)-(l27T) can 
be summarized in terms of the following algorithm. 

Given strictly feasible p = (pi,p2, ■■■,Pn)> t — £ » A* > 1, and e > 0, where t > is the 
step parameter of the barrier function method, fi is the step size of the algorithm, and e is the 
allowed duality gap (accuracy parameter) (see [fT8l ). do the following. 

1. Compute p(t) = argmin ||/ M g(9, p) d0-£ i=1 logft - J2 i= i log(p™ ax - Pi)j subject 
to Xli=i PiPi = PR starting at current p using Newton's method. 

2. Update p := p(i). 

3. Stopping criterion: quit if N/t < e 

4. Increase t — /it and go to step 1. 

In this algorithm, p(£) is called the central point and the first step of the algorithm is called the 
centering step. Then, at each iteration the central point p(i) is recomputed using Newton's method 
until N/t < e, that guarantees that the solution is found with accuracy e, i.e., P e (p)— P e (p op t) < e. 

For solving the optimization problem in the centering step, which is a convex problem with 
a linear equality constraint, the extended Newton's method is applied. It can be summarized as 
follows. 

Use p as an initial point and select error tolerance used in the Newton's method e r . 

1. Compute the Newton step Ap and decrement A(p). 

2. Stopping criterion: quit if A(p) 2 /2 < e r . 

3. Line Search: Choose step size s by backtracking line search. 

4. Update p := p + sAp and go to step 1. 

The Newton's step Ap and decrement A(p) are obtained from the following equations 

V 2 tf (p) a\ (*p\( -V# (p) 
a T U 
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where 



t 



I 



M-l 



h( P ) 



g(9,p) d6 



TV 



N N 



^Tiogp.-^iogbr*-^) 



and V 2 iJ(p) and ViJ(p) are the Hessian and gradient of H(p), respectively. Finally, decrement 
is defined as 



First, note that the Lagrange multiplier v' is non-negative, otherwise the equation (|46|) implies 
that 7^, % — 1, • • • , N are all positive numbers and, in turn, equation (1441) implies that pi = p™ ax , 
i = 1, . . . , N. Since it was assumed that Yli=i (3iPT ax > P-R> me condition (l48l) can not be satisfied 
for negative z/. As a result, z/ must be non-negative. Depending on whether bi is greater or smaller 
than v', two cases are possible. If bi < v' , then the condition (1451) holds true only if pi = 0. 
Indeed, if b t < v' , the expression fai/ + 7/ - (A(l + &iPt) + (1 - + friP;) 2 ) equals to 
the non-negative quantity + 7/ — when p; = 0. Furthermore, by considering the fact 
that the term fabi/ (A(l + + (1 ~ A)(l + hPi) 2 ) is strictly decreasing with respect to pi, 
it is resulted that the expression /3jz/ + 7/ — /^/(/^(l + 6^) + (1 — + hpi) 2 ) is greater 
than zero if p ri > 0, which means that the condition (|45T) can not be satisfied if p, > 0. Thus, Pi 
is zero when b{ < v' . Furthermore, if b{ > v' , then p, can not be equal to 0. It is because if pi 
equals zero then the condition (|44|l implies that 7/ = 0. Then, substituting Pi = and 7/ = in 
the condition (l46l) yields that 6, < z/, which contradicts the condition b{ > v' . Therefore, using 
the condition (|45T) and the fact that pi must be positive if bi > v' we can infer that 



A(p) = (Ap r V 2 tf ( P )Ap) 
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and a = (J3 U . . . 



Appendix C 
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By solving (l56l) with respect to p^, we can find that 
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Note that the negative root is not considered because Pi must be non-negative. We also defined 
above the function /j(z/,7-) for notation simplicity. By substituting the latter expression for p i 
into dUl) and (gT), we obtain 



If /j(z/,0) > pi max , then the function 7.) is strictly decreasing with respect to 7^ since v' 
is non-negative. Thus, the condition (1591 holds true only if 7^ > 0. In addition, the condition 
(1581) implies that p.; = pi max . Then the only remaining case is when < fi(v',0) < pi max . In 
this case, the conditions (l58l) and (l59l) hold true only if 7/ = 0. Considering the fact that the 
case of bi < v' is equivalent to the case when fi(v',0) < 0, we can conclude that all possible 
cases are analyzed and they then can be summarized as (|49l . □ 
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Fig. 2. Average SEP obtained through Monte-Carlo simulations and by using the closed form SEP expression l !24t . 
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Fig. 3. Average SEP obtained based on the proposed exact and approximate power allocation methods and the equal power 
allocation method with 3 relays. 
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Fig. 4. Average SEP obtained based on the proposed exact and approximate power allocation methods and the equal power 
allocation method with 5 relays. 
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Fig. 5. Average SEP obtained based on the proposed approximate power allocation method and the method of 1111 using 3 
and 5 relays. 
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